1/1 


SD-A143  268 
UNCLASSIFIED 


SIMPLE  EXAMPLES  OF  NON  LINEAR  TRUSS  8EHAVIOR(U> 
MINNESOTA  UNIV  MINNEAPOLIS  DEPT  OF  AEROSPACE 
ENGINEERING  AND  MECHANICS  P  G  HODGE  HUG  83  AEH-H2-2 
N88814-75-C-8177  F/G  28/11 


AD-A143  368 


\ 

w  < 


Report  AEM-H2-2 


SIMPLE  EXAMPLES  OF 
NON  LINEAR  TRUSS  BEHAVIOR 
by 

Philip  G.  Hodge,  Jr. 
Professor  of  Mechanics 


Department  of  Aerospace  Engineering  and  Mechanics 
University  of  Minnesota 
•  Minneapolis,  Minnesota  55455 


August,  1983 


Technical  Report 


j  I  fltwAi&iad  w>ay— otefett**  ''np-*  *■■■' 

i 

t 

Prepared  for 

OFFICE  OF  NAVAL  RESEARCH 
_  (  Arlington,  VA  22217 


OFFICE  OF  NAVAL  RESEARCH 
Chicago  Branch  Office 
536  South  Clark  St. 
Chicago,  IL  60605 


l!»i3  <t^  uijunal  k.ua  inti 
to  puKtto  r»i«uKi  .rod  igUj  If 

dtoMbartaa  it  mMmlrnd 

' —  — .  ,i 


(jm 


nnn 


84  07  12  009 


.  _ •  _ . .  ^  * 

SECURITY  CuA^SI^ICATioh  Of  ThiI  9* ACC  B*l***<> 


REPORT  CCCUMEHTATICH  PAGE 


i  hipo^t nuuttn 

AEM-H 2-1- 

4.  TIT Lt  f»4  JOWH.J 


S*rc  READ  D43TX’JCT10NS 

_ beei-.re  complettng  rorm 

;j.  ooyr  accession  noJ  >•  recipient's  catalog  hTTmTT^ 


SIMPLE  EXAMPLES  OF  NON  LINEAR  TRUSS 
BEHAVIOR 


••  TYPE  O  F  4l»O^T  4  pjwoo  COVEHEO 

Technical  Report 

8.  0*a.  4!PO»T  NUU«!P 


7.  authors 

P.G.  HODGE,  JR. 


-ohT^a^T  oa  orant  wuCaeSTIJ 

M14-75-C-0177 


*  PERFONMINO  ORGANIZATION  NAME  ANO  4DORI13 

University  of  Minnesota 
Minneapolis,  Minnesota  55455 


PAOO  PAM  El.  !«WT,  PPOJECT  TA3< 
ARIA  A  WORK  UNIT  HUMIC  AS 


MF  064-429 


I II.  COMTROLLINO  OFFICE  NAME  AMD  ABDREJS  12.  REPORT  OATS 

OFFICE  OF  NAVAL  RESEARCH  ’August.-,  198  3 _ 

Arlington,  VA  22217  '*•  Humm€n  PMa* 

_ _ _ 29  _ 

1 4.  monitoring  agency  NAME  8  AOORXSV"  IIMmwI  trom  CanlrmlllKg  Olll  n)  13.  EEC  UNITY  CLASS.  JZt  tfilm 


OFFICE  OF  NAVAL  RESEARCH 
Arlington,  VA  22217 


OFFICE  OF  NAVAL  RESEARCH 
Chicago  Branch  Office 
536  South  Clark  St. 

Chicago.  II.  606  OS  . 

18.  01JTRI9CTION  STATEMENT  (ot  UUmiXfQup  -  x  .  ..  u. ...  ,.^a  C$>’r*luv 

iot  psiOc  n.i.-tm  >-?o£\  solo:  Its 

I  dkrtrttu^ou  U  us&aitted. 


Unclassified 

(»«.  OZCLASSiFICATION/DOOMONAOInO 
tCMeDULE 


17.  OISTNISUTION  JT ATEMSM T  (•»  S  «  •IiIihI  Mlml  I*  Slot*  39,  It  dltlortnl  trrm  Rfyort) 


IS-  SUPPLEMENTARY  NOTES 


•  9-  KEY  WOAOJ  fC«n<itNM  on  r*v*r««  »•<#•  II  wocMiay  m>4  I4m%tlfy  tty  WocJl  pylwj 


“Plasticity,  buckling,  nonuhiquehess t  unloading,  superpositic 


20.  y  ACT  (Ctmtlmf  cn  •*  f  ll  r**o»*«**7  **4  Hmrttfy  by 

"Nonlinear  material  models  are  defined  to  represent  elastic/ 
perfectlv-plastic  and  elastic/buckling  behavior.  A  simple  three- 
bar  truss  is  used  to  demonstrate  that  under  a  monotonical ly 
increasing  prescribed  displacement  the  truss  nav  exhibit  reverse 
stressing  or  non-uniqueness,  and  that  when  two  different 
control  displacements  are  applied  the  principal  of  superposition 
does  not  hold.  , 


00  I  j°m“,  1473  Winaii  Of  I  »OY  11 11  o»Mi*T8 

</»  010J-0U-#»01  I 


88CURITY  Clarification  of  thi*  rasa 


I 


SIMPLE  EXAMPLES  OF 


NONLINEAR  TRUSS  BEHAVIOR 


Philip  G.  Hodge,  Jr. 


Professor  of  Mechanic: 


University  of  Minnesota 


1.  Introduction.  The  theory  of  linear  elasticity  is  a  very 
"nice"  mathematical  theory  with  many  convenient  features  such  as 
superposition,  uniqueness,  and  the  equivalence  of  proportional 
loading  to  proportional  stressing.  Theories  involving  nonlinear 
material  behavior  do  not  necessarily  have  these  nice  features,  and 
it  is  instructive  to  consider  very  simple  examples  which  dramatically 
illustrate  this  fact. 

In  1951  Drucker  (1)  considered  a  simple  3-bar  truss 
essentially  similar  to  the  one  shown  in  Fig.  1  which  is  subjected 
to  a  single  monotonically  increasing  vertical  load  corresponding 
to  the  load  Q.  In  Drucker' s  example  each  bar  is  made  of  an 
elastic/perfectly-plastic  material  with  identical  elastic  properties 
but  different  yield  strengths.  For  a  suitable  choice  of  yield 
strengths,  bar  1  will  first  yield  in  compression  but  will  then 
unload  and  will  eventually  yield  in  tension  when  the  truss  fails.  ' 
The  present  author  (2,  3)  has  made  frequent  use  of  trusses  similar 
to  the  one  in  Fig.  1. 

The  present  note  shows  how  this  truss  can  be  used  to 
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discuss  two  different  models  of  non-linear  material  behavior 
in  relation  to  the  above-mentioned  features  of  superposition, 
uniqueness,  and  proportional  stressing. 

The  two  models  to  be  treated  are  the  familiar  elastic/ 
per  feet  ly -plast  ic  (K/PP)  one  illustrated  by  t  he  dashed  curve-  in 
Fig.  2,  and  an  idealized  elastic  buckling  (K/b)  model  indicated 
by  the  solid  curve.  Since  we  are  primarily  concerned  with 
compressive  stresses,  we  define  the  bar  shortening  and  the  neg¬ 
ative  bar  force  by 
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(1) 


The  buckling  model  places  no  limit  on  negative  values  of 

* 

(it  would  be  a  trivial  extension  to  construct  a  combined 
model  wnich  yielded  in  tension  and  buckled  in  compression)  but 
an  equally  important  difference  in  the  two  models  occurs  on 
unloading.  At  point  B,  for  example,  the  t/PP  model  unloads 
along  BCD  whereas  the  E/b  model  retraces  the  loading  path  baog. 

The  next  section  lists  the  defining  equations  for  the  truss 
in  Fig.  1  according  to  the  two  models,  and  the  concluding  section 
defines  three  specific  examples  chosen  to  illustrate  reverse  stressing 
under  a  monotonic  load,  non-uniqueness,  and  superposition,  in 
that  order.  A  summary  of  results  are  given  in  Section  3,  and 
some  details  of  the  solutions  are  included  in  the  Appendix. 

2.  Equations.  The  cross  bai  AB  is  assumed  rigid  and  the 
three  vertical  bars  all  have  the  same  modulus  E,  area  A,  and 
minimum  moment  of  inertia  I.  We  find  it  convenient  to  use 
asterisks  to  represent  physical  quantities  and  define  dimensionless 
variables  as  indicated  below.  For  both  models  we  define  the  still 


ness  of  bar  2  by 


1 


k  =  AE/H  ( 2 ) 

where  H  is  the  length  of  the  bar.  For  the  E/lt  model  we  use  the 
simple  Euler  buckling  formula  for  a  pinned  bar  and  denote  t  lie 
shortening  of  bar  2  at  the  onset  of  buckling  by 
C  2rT  2, 

-  crnEIH_nI  ,  , , 

3  =  "k - “ T  ’  AE  "  AH_  (J) 

We  then  define  dimensionless  displacements  v  and  w,  loads  1’  and  (j, 
shortenings  s^, and  compressive  forces  by 


*  __ 

V  =  V  /s 

#  _ 

w  -  w  /s 

*  _ 

s  =  s . /s 

1  1 

*  _ 

*  _ 

*  — 

P  =  P  /ks 

Q  =  Q  /ks 

C  =  C . /ks 

l  i 

For  the  E/PP  model,  let  V  be  the  yield  force  in  bar  2  and 
define  the  above  variables  by  the  similar  formulas 

•  •  * 

v  =  kv  /Y  w  =  kw  /Y  Bi  =  ksi/Y 

(5) 

P  =  P*/Y  Q  *  Q*/Y  Ci  =  C*/Y 

We  denote  the  length  of  bar  i  by  cmH  where  u2  =  1  and  we 

will  choose  and  in  each  example.  In  view  of  the  above 
definitions,  the  shortening  at  whi  h  bar  i  buckles  is  then  given 
by 

si  =  l/ui  (>>) 
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So  that  the  two  models  will  have  the  same  scale  in  Fig.  2,  we 

» 

choose  the  yield  stress  in  bar  i  to  be  inversely  proportional  to 
the  square  of  its  length,  hence 

Yi  -  1/a.2  (7) 

We  assume  that  all  displacements  are  small  so  that  the  only 
non-linearity  is  the  material  behavior,  and  the  three  bars  can  be 
taken  to  remain  vertical.  The  kinematics  are  given  by 

6si  =  (4Aw  -  3Av,  Av,  2Av  -  Aw)  (8) 

where  we  find  it  convenient  to  use  the  incremental  form  for 
all  equations.  The  statics  are  obtained  from  moment  equilibrium 
about  the  two  loads: 

AP  =  -3AC1  +  AC2  t  2AC3  (9a) 

AQ  =  4AC1  -  AC3  (9b) 

Finally,  the  constitutive  equation  for  each  bar  in  the  E/PP 
model  may  be  written 

IF  (Ci  ~  ♦l/ai2  AND  CiAsi  >  0) 

THEN  AC3  =  0  ELSE  AC^  =  ASj/Oj  (10) 

For  the  E/B  model  the  equation  is  similar  but  significantly 
different: 


IF  s.  >  1/a.  THEN  AC.  =  0 
l  l  i 


=  AS./ui 


( 


i 


( 


( 


ELSE  ACV 


(11  ) 


b 

In  each  case  the  ELSE  clause  represents  elastic  behavior  and 
the  IF  clause  represents  the  Idealized  inelastic  behavior.  When 
all  three  bars  are  elastic,  the  solution  for  either  mode]  may  be 
written 


AP 


AQ 


(  —  +  1  +  l_)Av 
al  a3 


(ii  ♦  —  - )  Aw 
al  a3 


-(  —  +  — )Av  +  (—  +  i-)Aw 
al  ft3  al  a  3 


(12a) 


(12b) 


Equations  (12)  could  be  solved  for  Av  and  Aw.  However,  in  each 
of  our  examples  we  will  apply  only  one  load  at  a  time  so  that  the 
zero  load  equation  is  trivially  solved  to  relate  Av  and  Aw  and  the 
other  equation  relates  load  and  displacement.  In  all  cases, 
we  will  regard  one  of  the  displacements  as  the  control  variable. 

3^ « _ Examples .  As  an  example  to  illustrate  unstressing 

we  take  =  (2,  1,  4),  let  Q  be  the  only  load,  and  increase  the 
control  displacement  w  under  Q  from  0  to  3.  For  the  E/PP  model 
this  example  is  essentially  similar  to  Drucker’s  [1). 

The  results  are  shown  in  Table  1  and  Figs.  3  and  4;  details 
are  given  in  the  Appendix.  The  "status"  column  in  Table  1 
shows  for  each  bar  if  it  is  elastic  (E) ,  yielding  in  tension  (T) 
or  compression  (C)  ,  or  buckling  (B) .  A  "stage"  is  the  time  spent 
with  no  bar  changing  status,  stage  1L  is  the  limit  of  stage  1  as 
bar  1  reaches  yield  and  changes  from  E  to  C,  etc. 

As  shown  in  Table  1  as  w  is  increased  first  bar  3  and  then 


bar  1  yield  in  compression.  However,  as  is  clear  from  Fig.  1, 
a  mechanism  motion  with  bar  2  rigid  would  require  bar  3  to 


I 


lengthen.  Therefore,  in  stage  3  bar  1  elastically  unloads  through  zero 
and  eventually  yields  in  tension  as  the  yield-point  load  is  reached. 

The  dashed  curves  in  Pigs.  3  and  4  respectively,  show  the  load- 
displacement  history  and  the  shortening  history  of  bar  3. 

Stages  1  and  2  do  not  involve  unstressing  or  tensile  yield  so 
that  they  are  the  same  tor  the  L/B  model.  llowevei  ,  ir;  stage  3 
as  bar  3  starts  to  reclaim  its  buckl  ing  deformation,  there  is  no 
change  in  its  force.  Since  bar  1  is  now  buckling  also  re¬ 
mains  constant.  Therefore,  equilibrium  shows  that  C  and  the  load 
are  also  constant.  Thus,  as  w  increases  bar  AB  rotates  about  the 
unchanged  location  of  the  top  of  bar  2.  However,  this  situation 
lasts  only  until  Sj  is  reduced  to  at  which  point  it  resumes 
elastic  behavior  in  stage  4.  Since  tensile  yielding  is  not 
considered,  stage  4  continues  until  bar  2  reaches  buckling.  The 
buckling  collapse  mechanism  of  the  truss  thus  involves  a  rotation 
about  the  top  of  bar  1.  The  results  are  shown  by  the  solid 
curves  in  Pigs.  3  and  4.  Notice  the  finite  horizontal  portion  of 
the  load-displacement  curve  in  stage  3,  followed  by  further  increase 
of  load  in  stage  4. 

* 

As  a  second  example  to  illustrate  non-uniqueness,  we  consider 
a  truss  with  ,r  -  (4,  1,  2)  where  P  is  the  only  load.  The 
control  variable  v  is  to  be  increased  from  0  to  1 .  Pot  the  E/pp 
model  the  all-elastic  stage  1,  shown  in  the  top  line  of  Table  2 
ends  when  bars  1  and  3  both  reach  compressive  yield  at  the  same 
instant.  Therefore,  in  stage  2,  AC^  -  AC^  =  0  and  since  there  is 
no  load  U  Kq .  (9b)  becomes  an  identity.  Therefore,  tire  only 

information  about  Aw  comes  from  (8)  and  the  inequalities  in  (ID) 
f or  bai s  l  and  3 : 

•Other  simple  examples  of  non-uniqueness  ate  discussed  in  14)  and  )  ‘>  |  . 
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As^  =  4 Aw  -  3 Av  >  0  As^  =  2Av  -  Aw  >  0  (ii) 

These  must  apply  for  any  infinitesimal  inclement  in  stage  2, 
which  leads  to 

3/4  <  dw/dv  <2  (14) 

as  shown  in  the  last  column  for  stage  2.  During  this  stage  the 
load  P  must  increase  with  v  and  the  bar  forces  are  all  unique, 
but  w  may  take  any  value  permitted  by  (14).  However,  we  note 
that  once  w  has  been  established  for  any  particular  v,  the 
restrictions  (14)  apply  from  that  point.  Therefore,  althoug! 
not  unique  as  v  is  increased,  any  solution  reached  is  stable 
if  v  is  held  constant  at  any  time. 

Figure  5  shows  the  strain-path  trajectory  of  bars  1  and  3. 
rt  is  unique  along  OA  in  stage  1,  but  during  stage  2  it  may  follow  any 
path  with  positive  slope  in  the  domain  ABC.  In  particular, 
at  stage  2L  it  may  have  reached  any  point  on  the  line  BC. 

However,  if  the  solution  at  v  =  3/4  is  observed  to  be  at  point  D,  say, 
then  the  possible  solutions  when  v  =  1  are  restricted  to  the  segment 
F.F . 

For  the  E/B  model  the  equations  and  unique  part  of  the 
results  are  exactly  the  same,  so  that  we  have  not  repeated  them 
in  Table  2.  However,  the  inequalities  apply  to  the  total  shortening 
rather  than  instantaneous  increments,  so  that  (13)  and  (14)  must 
be  replaced  by 

sL  >  s1  =  1/4  s2  >  s2  =  1/2 

(3/4)v  +  1/16  <  w  <  2v  -  1/2  (IS) 


ft 


« 


1 


« 


A 
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Not  only  is  the  solution  for  w  not  unique,  it  is  only  neutrally 
stable  and  could  change  from  one  value  to  another  with  no  change  in 
v.  Thus,  at  stage  2L  it  could  be  any  point  on  BC  regardless  of 
its  earlier  values  at  v  =  3/4. 

The  final  example  discusses  superposition.  We  return  to 
the  first  truss  where  =  (2,  1,4)  and  control  both  v  and  w, 
increasing  them  from  zero  to  w  =  1/4,  v  =  1.  For  the  E/PP 
model  the  final  loads  depend  upon  the  order  in  which  the  displace¬ 
ments  are  increased.  Suppose  that  they  are  applied  in  the  order 
v,  w;  i.e.,  v  is  increased  to  1  with  w  held  at  zero,  then  v  is 
held  at  1  while  w  is  increased  to  1/4.  The  light  solid 
curve  ABCDE  shows  the  history  of  load  Q;  P  would  have  a  similar 
curve.  Complete  results  may  be  found  in  the  Appendix.  The 
final  state  when  v  =  1,  w  =  1/4  is  shown  in  line  1  of  Table  3. 

On  the  other  hand,  if  w  is  first  increased  to  1/4  and  then  v 
is  increased  to  1,  the  history  of  Q  is  given  by  the  heavy  solid 
curve  AFGI1IJ  in  Fig.  6.  Not  only  is  the  history  quite  different, 
but  the  final  load  values  (points  E  and  J)  do  not  even  agree  in  sign. 
Line  2  of  Table  3  lists  all  final  values. 

If  the  principle  of  superposition  were  valid,  we  could 
find  two  components  corresponding  to  v-only  and  w-only  displace¬ 
ments  and  then  add  them.  Lines  3  -  5  in  Table  3  show  the 
final  results,  and  we  see  that  line  5  is  quite  different  than  either 
lines  l  or  2.  In  terms  of  Fig.  6,  we  could  add  the  w-only  “elu¬ 
tion  AFG  to  the  v-only  solution  ABCD  by  translating  the  first 
curve  so  that  A  is  at  D.  The  resulting  curve  ABC ON M  has  its 
new  part  shown  light  dashed. 

The  order  of  superposition,  of  course,  does  not  metier. 

It  v-only  is  added  to  w-only,  the  curve  AFGKI.M  has  exactly  the  s  mil- 
terminal  point  M. 


For  the  E/B  model  the  shortenings,  bar  forces,  and  loads  are 
all  unique  functions  of  the  instantaneous  displacements,  so  that 
curves  ABCD  and  AFGHD  in  Fig.  7  both  end  at  the  same  state  0. 

The  complete  solution  at  the  final  point  v  =  1,  w  =  1/4  is  shown 
in  line  6  of  Table  3.  Observe  that  the  loads  are  very  much 
different  from  those  required  by  the  E/PP  model,  essentially 
because  the  latter  had  bar  1  with  substantial  yielding  in  tension. 

However,  as  shown  by  lines  7,  8,  9  in  Table  3  and  by 
curves  ABCL  or  AFKL  in  F^g.  7,  the  results  of  superposition 
still  do  not  agree  with  the  actual  solution. 


10 


APPENDIX 


We  present  here  some  of  the  intermediate  steps  ul  t ) j •  • 
examples  in  Section  3.  Durincj  a  staye  in  which  all  burs  aie 
elastic,  it  follows  from  Eqs .  (8)  and  the  "ELSE"  part  of  (lu) 
or  (11)  that 


AC. 

l 


,4 Aw  -  3Av 


Av, 


2AV  -  Aw. 


1 A  1  ) 


Further,  if  any  bar  is  inelastic,  the  corresponding  component  of 
AC^  is  replaced  by  zero.  Therefore,  we  can  immediately  write 
an  explicit  expression  for  AC^  in  terms  of  Aw  and  Av  foi  any 
staye  of  loadiny. 

Example  1.  In  this  example 


«i  =  (2,  1,  4)  Yt  =  (4,  16,  1 ) / 1 6 

so  that  Eq .  (Al)  can  be  written 


(2,  4, 


( A  2  > 


AC\  =  (8,  0,  - 1 ) Aw/ 4  +  (-3,  2,  1 ) Av/2  (A3) 


The  only  non-zero  load  is  Q,  hence  it  follows  from  Eq.  (9a) 

AP  =  -3AC1  +  AC2  +  2AC’3  =  0  (A4) 

Substitution  of  (A3)  or  its  partially  inelastic  t epl aoemont  in 
( A4 )  produces  an  equation  which  is  easily  solved  for  Av  in  terms 
of  the  control  variable  Aw.  Thus  we  beqin  with 
St aye  K  EKE 

AP  =  (-24-2)  Aw/4  ♦  (9  +  2  +  2)  Av/2  =  () 


1 ) 

For  the  E/PP  model  this  stage  will  end  when  bar  3  reaches 
compressive  yield,  i.e.,  when 

C3  =  0  +  Aw/4  =  1/16  Aw  =  w  =  1/4  (A6) 

Equations  (A5)  ,  (A6),  and  (8)  -  (11)  then  determine  the  complete 
solution  fox  stage  II,  as  given  in  line  1  of  Table  1. 

For  stage  2,  Eq.  <A5)  is  replaced  by 
Stage  J2  .  EEC 

AC.  =  (2,  0,  0) Aw  +  (-3,  2,  C)Av/2 

AP  =  -6 Aw  +  (9  +  2 ) Av/2  =0  Av  =  (12/11) Aw  (A7) 

This  stage  ends  when  bar  1  yields  in  compression.  Since  C^  =  1/8 
at  the  end  of  stage  1  (see  line  1  of  Table  1) ,  we  have 

C3  =  (1/8)  +  (4/11) Aw  *  1/4  Aw  =  11/32  (A8) 

Equations  (A7),  (A8) ,  and  (8)  -  (11)  then  give  the  complete 

increment  solution  during  stage  2: 

Av  =  3/8  As i  =  (1/32) (8,  12,  13) 

ACa  =  (1/32) (4,  12,  0)  AQ  =  1/2  (A9) 

When  these  are  added  to  the  stage  1L  values  in  line  1  of  Table  1, 
we  obtain  the  complete  solution  at  stage  2L  as  shown  in  line  2. 

At  this  point  it  might  appear  logical  to  assume 
Stage  3X.  CEC 

AC3  =  (0,  1,  0) Av 


AP  =  Av  =  0 


Asa  =  (4,  0,  - 1 ) Aw 


( A 1  0  > 
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However,  since  Aw  is  positive,  we  would  be  predicting  <  o 

which  is  not  consistent  with  the  assumption  of  yielding.  Thoie- 
fore,  instead  we  write 
Stage  2-  CUE 

ACa  =  (0,  0,  -1 ) Aw/ 4  +  (0,  2,  l)Av/2 

AP  =  -Aw/2  +  (2+2)  Av/2  =  0  Av  =  Aw/4  (All) 

This  stage  ends  when  bar  3  yields  in  tension : 

C3  =  1/16  -  Aw/8  =  -1/16  ' 

Aw  =  1  Av  =  1/4  4  s.  =  (13,  1,  -2)/4 

ACt  =  (0,  2,  -l)/8  AQ  =  1/8  (A12) 

Addition  of  these  values  to  those  in  line  2  of  Table  1  produces 

the  values  in  line  3. 

in  the  final  stage 
Stage  4 .  CET 

CL  =  (0,  Av,  0) 

AP  =  Av  =  C^Ci  =  0  Ast  =  Aw  ( 4 ,  0,  -1)  (M3) 

which  is  the  collapse  mechanism  about  the  end  of  bar  2.  Since 
our  program  calls  for  w  to  increase  to  3,  we  set  Aw  =  45/32 
in  Eqs .  (A13)  and  add  to  line  3  to  get  line  4  in  Table  1. 

For  the  E/B  model,  Eqs.  (A5)  througti  (A9)  are  still  applicable 
hence  the  solution  is  exactly  the  same  through  stage  2h. 


1  3 

However,  since  s^  =  21/32  is  beyond  the  buckling  value  =  1/4, 

Eq.  (A10)  is  a  valid  description  with  bar  3  unstressing,  but 
remaining  buckled,  until 

s j  =  21/32  -  Aw  =  1/4  Aw  =  13/32 

Av  =  AP  =  AQ  =  ACi  =  0  As.  =  (13/32) (4,  0,  -1)  (A14) 

Addition  of  these  values  to  line  2  gives  line  5  in  Table  1. 

Bar  3  now  resumes  elastic  behavior,  hence  stage  4  is  described 
by  Eq.  (A8) .  Since  tensile  yield  is  not  considered,  this  stage 
continues  until  bar  2  buckles: 


s2  ^  5/8  +  Aw/4  =  1 

Aw  =  5/2 

Av  =  5/8 

ACa  =  1/16(0,  6,  -3) 

AQ  =  3/16 

As  ^  =  (1/8)  (39,  3,  -6) 

(A15) 

Addition  of  (A15)  and  line  5  of  Table  1  produces  line  6. 

In  the  final  stage 
Stage  5.  BBE 

AC\  =  (0,  0,  2 Av  -  Aw)/4 

AP  =  Av  -  Aw/2  =0  Av  =  Aw/2 

AQ  =  6Ci  =0  ASi  =  (5,  1,  0)Aw/2  (Alb) 

This  is  a  mechanism  motion  of  rotation  about  bar  3  and  will 
continue  until  w  reaches  its  final  value  of  Thus  line  7  ut 

Table  1  is  obtained  by  setting  Aw  =  1/2  in  (Alb)  and  adding 
the  result  to  line  6. 
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Example  2.  In  this  example  Eqs.  (A2)  -  (A4)  are  replaced  by 


«i  =  (4.  1,  2) 


Y  =  (1,  16,  4)/lb 


si  -  (1.  4,  2 ) / 4 


AU  -  4AC^  -  AC  i  ■-  0 


( A 1  V  ) 

(Alb) 


AC i  =  (2,  0,  - 1 ) Aw/2  +  (-3,  4,  4)Av/4 


(Alb) 


Therefore,  we  can  find  Aw  in  terms  of  the  control  variable  Av. 
When  all  bats  are  elastic,  we  use  (A19)  and  obtain 
Stage  1.  EEE 

AQ  =  (8  +  1 ) Aw/ 2  +  (-12  -  4)Av/4  =  0 


Aw  =  (8/9) Av 


(A20) 


For  the  E/PP  model  this  stage  ends  when  bars  1  and  i  reach  yield 
simultaneously  at  v  =  Av  =  9/20  which  leads  to  the  values  in 
line  1  of  Table  2. 

In  stage  2  bars  1  and  3  are  both  yielding  hence 
Stage  2.  CEC 


AC^  =  (0,  1,  0) Av 


AQ  -  0  AP  -  Av 


Asl  =  Aw (4 ,  0,  -1)  +  Av ( - 3 ,  1,  2) 


(A21 ) 


The  only  information  available  for  Aw  is  in  the  Inequa l  i t ies 
(13).  This  stage  ends  when  the  remaining  bar  yields: 


1  5 


C2  =  9/20  +  Av  =  1  Av  =  11/20  (A22) 

With  this  value  of  Av,  the  sum  of  liues  1  and  2  in  Tabic  i 
g i ves  line  3  . 

Example  3.  In  these  final  examples  the  truss  properties 
are  again  given  by  (A2)  and  the  fully  elastic  force-increment 
solution  by  (A3).  Both  v  and  w  are  controlled,  and  we  will  find 
it  convenient  to  express  the  kinematic  and  static  equations 
(8)  and  (9)  in  integrated  form: 

sx  =  (4,  0,  —  1 ) w  +  (-3,  1,  2 ) v  ( A 2 3 ) 

0  -  «cA  -  c3 

P  =  -3C1  +  C2  +  2C3  (A24) 

Thus,  given  the  control  variables,  (A23)  gives  an  explicit 
expression  for  the  shortenings.  Using  (A2)  or  its  non-elastic 
replacement,  we  find  the  force  increments,  but  we  can  wait 
and  use  the  integrated  forces  to  find  the  loads  directly  from  (A24). 

In  example  3A  we  first  increase  v  to  1  and  then  increase 
w  to  1/4;  in  3B  we  reach  the  same  final  values  by  increasing 
first  w  and  then  v. 

The  computations  fall  into  a  simple  pattern  and  are  conveniently 
presented  in  tabular  form  in  Table  4.  We  consider  first  example 
3A  tor  an  E/PP  material.  The  zeroes  in  the  unnumbered  top 
line  emphasize  that  we  start  from  the  zero  state  and  obtain  the 
solution  by  finite  increments. 
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In  stage  1  all  bars  are  elastic  and  v  is  the  active  coni  i  ol 
variable.  We  denote  its  increment  from  zero  by  A.  in  line  1, 
column  1,  and  express  the  bar-force  inciements  in  teuus  of 
in  column  3.  In  this  incremental  stage  we  do  not  need  the  i terns 
in  columns  2  ot  4.  Stage  1  will  terminate  when  bar  i  i  caches  i  i 
yield  value  in  compression.  The  resulting  equation  for  the 
increment  A  is  written  in  column  5,  and  its  solution  is  enteiod 
in  column  6. 

The  solution  at  stage  1L  is  then  written  in  line  2. 

Values  for  the  displacements  in  column  1  and  bar  forces  in  column 
3  are  obtained  by  putting  the  value  A  =  1/8  in  line  1  and 
adding  the  result  to  the  zero  values  at  the  beginning  of  stage  1 . 
Although  we  could  treat  the  shortenings  s  and  loads  Q  and  p 
in  the  same  incremental  fashion,  it  is  more  convenient  to  use 
( A 2 3 )  and  (A24)  and  complete  columns  2  and  4  directly  from  the 
w,  v,  and  Ch  values  in  line  2. 

Line  3  treats  the  next  increment  in  the  same  way.  In  column 
3,  AC^  and  aC2  still  have  their  elastic  values  from  (A3),  but 
AC^  =  0  since  it  is  already  yielding.  Stage  2  terminates  when 
bar  1  yieids  in  tension,  which  leads  to  the  equation  in  column  5 
and  the  solution  A  =  1/24  in  column  6. 

For  line  4  we  set  A  =  1/24  in  line  3  and  add  the  result,  to 
line  2  to  obtain  values  of  w,  v,  and  Ch  .  Again,  the  values  of 
s^,  Q,  and  P  are  obtained  directly  from  (A23)  and  (A24). 

Lines  5  and  6  contain  the  solution  for  stages  3  and  31. 
figuied  in  the  same  way.  No  bar  changes  state  before  v  reaches 
its  final  value,  hence  the  equation  in  column  5  is  simply  v  1. 
Line  6  is  also  recorded  as  line  3  m  Table  3  as  the  v-component 
for  "superposition''. 
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In  stage  4  control  switches  to  w.  Since  the  elastic- 
solution  for  w  increasing  produces  a  positive  AC^ 

and  negative  both  plastic  bars  now  unload  and  all  three  bars 

are  elastic  in  line  7.  Now  w  reaches  its  final  value  with  no  bar 
changing  state  so  the  equation  in  column  5  is  w  =  1/4.  Line 
8  gives  the  final  values  for  example  3A;  the  same  values  are 
recorded  in  line  1  of  Table  3. 

Example  3B  for  the  E/PP  material  is  treated  in  exactly 
the  same  way.  The  details  are  displayed  in  lines  9-18  of  Table  4, 
and  the  final  results  from  line  18  are  recorded  as  line  2  of 
Table  3.  Also,  line  12  of  Table  4  is  repeated  as  line  4  of 
Table  3  for  the  w-component  for  "superposition". 

For  the  E/B  model  and  example  3A,  stages  1,  1L,  and  2  are 
the  same  as  for  the  E/PP  model.  However,  since  the  E/B  model 
does  not  have  any  constraint  on  allowable  tensile  loads,  the 
terminating  condition  is 

v  =  1/8  +  A  =  1  A  =  7/8  (A25) 

Substitution  of  this  value  in  line  3  of  Table  4  and  addition 
to  line  2  gives  the  stage  2L  solution  in  line  19  for  the  E/B 
model.  This  line  is  also  written  as  line  7  of  Table  3. 

Since  bar  3  has  undergone  a  finite  amount  of  buckling, 
it  will  remain  in  the  buckled  state  in  stage  3  as  shown  in  line 
20  of  Table  4.  In  fact,  it  is  still  buckled  when  w  reaches  its 
final  value  of  1/4.  The  final  values  at  stage  31.  are  listed  in 
line  21  of  Table  4  and  line  6  of  Table  3.  From  column  2  and 
the  last  Eq.  (A2)  we  verify  that  bar  3  is  still  well  within  the 
buckled  state,  bar  1  is  in  tension,  and  bar  2  is  on  the  verge  ol 


buckling. 
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In  example  3B  the  two  models  are  exactly  the  same  through  stage 
2L.  In  stage  1  bar  1  will  remain  in  the  buckled  state  until 
its  shortening  reduces  to  .  To  determine  buckling  or  unbuckling 

criteria  in  this  example  it  is  necessary  to  have  the  As^ 
expressions  available  as  given  for  example  in  line  22  column  2. 
in  this  case  the  stage  ends  when  bar  1  unbuckles  at  •  -  1/b. 

Line  23  shows  the  solution  at  stage  3L.  The  shortening  in 
column  2  is,  of  course,  the  same  whether  it  is  computed  by 
adding  the  increment  in  line  22  to  the  value  at  stage  2L  in  line  12 
or  directly  from  Eq.  (A23) . 

The  rest  of  the  table  is  completed  in  the  same  fashion.  As 
noted,  lines  27  and  21  are  identical  for  this  model. 


i 
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